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Abstract. We review how the use of recent precise data on kaon decays together with forward 
dispersion relations (FDR) and Roy's equations allow us to determine the sigma resonance pole 
position very precisely, by using only experimental input. In addition, we present preliminary 
results for a modified set of Roy-like equations with only one subtraction, that show a remarkable 
improvement in the precision around the a region. For practical applications, these results are shown 
to be very well approximated by a very simple conformal expansion. 
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INTRODUCTION 

The values quoted in the Particle Data Table for the sigma mass and width, based on 
both pole position and Breit-Wigner parameter determinations are very widely spread, 
with an estimated mass and half width of [1]: 

^ = Ma- iTa/2 ~ (400 - 1200) - /(250 - 500) (MeV). (1) 

This large uncertainty is mainly due to the fact that old data sets for pion-pion scattering 
are poor and often contradictory. Moreover, the choice of data sets varies among different 
works. To make things worse, there is quite a variety of different ways to extrapolate the 
data on the real axis to the complex plane, and the pole position of the sigma is greatly 
affected by model dependences. 

This said, model independent techniques for extrapolating amplitudes from the real 
axis onto the complex plane exist in the form of dispersion relations, which allow us 
to analytically continue an amplitude away from the real axis provided we know its 
imaginary part for physical values of the energy. These dispersive techniques have 
already been successfully used for predicting the position of the sigma pole, with a 
remarkable agreement among the different works: 

440 - HAS MeV Dobado, Pelaez (1997) lH (2) 

470 ± 50 - / 260 ± 25 MeV Zhou et al. (2005) [H (3) 
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In particular, there exists a dispersive representation that incorporates crossing exactly, 
written by Roy [4], which involves only the partial wave amplitudes. Roy's equations 
have already been used to predict the position of the sigma pole from the theoretical 
predictions of ChPT [5], obtaining: 

v^ = 441^16 _/272l?^ 5 (4) 

In addition, the data coming from the E865 collaboration at Brookhaven [@], and 
especially the recently published data from NA48/2 ff] provide us with very precise 
data on pion-pion scattering at very low energies. These allow us to obtain very reliable 
parametrizations of the SO wave at low energy from which the scattering lengths 
can be directly extracted [9] with a remarkable precision and in good agreement with 
the theoretical predictions of ChPT [5]. 

Our aim is thus to perform a dispersive analysis, including all available experimental 
data, in order to give a precise and model independent determination of the sigma pole 
position, by using exclusively data, analyticity and crossing symmetry. We use both 
Forward Dispersion Relations (FDR) and Roy's equations, without assuming ChPT, so 
that we can actually test its predictions. 



APPROACH AND RESULTS 

The details on the parametrizations used for the data have been explained fully in this 
conference (see talk by J. R. Pelaez (^), as well as in Ref. [10] - that we will denote by 
KPY08. It is enough to say here that two different sets of parameters are considered: 

• Unconstrained Fits to Data (UFD), in which each partial wave is fitted indepen- 
dently. This set satisfies both FDR and Roy's equations within the experimental 
errors in all waves except the Roy equation for the S2 wave, for which the devi- 
ation is about 1.3 a, and the antisymmetric FDR above 930 MeV by a couple of 
standard deviations. 

• Constrained Fits to Data (CFD), obtained by constraining the fits to satisfy simul- 
taneously FDR and Roy's equations, so that all waves are correlated. The CFD set 
provides a remarkably precise and reliable description of the experimental data, and 
at the same time satisfy the analytic properties remarkably well. 

These two sets provide a reliable parametrization for the imaginary part of the partial 
waves that we need as input for Roy's equations. 

An elastic resonance has an associated pole on the second Riemann sheet of the 
complex plane S-matrix, which, as it is well known, corresponds by unitarity to a zero on 
the first sheet. As usual then, we just need to look numerically for zeroes of the S-matrix 
on the physical sheet: 

S^o{s) = l+2ia{s)t^{s), (5) 



where the analytic extension of the partial wave amplitudes away from the real axis is 
given by Roy's equations: 
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The domain of validity of Roy's equations has been shown to cover the region of the 
complex plane where the sigma lies 

Taking the UFD set as the input for Roy's equations, we find an S-matrix zero at 
a/v = (426 ± 25) — ?'(241 ± 17) MeV. However, Roy's equations are not completely 
satisfied by this data set, thus the pole position will be much more reliable if the input 
satisfies the equations, as it is the case for the CFD set. In this case we find: 



which still has big uncertainties due to the strong dependence of Roy's equations on 
the scattering lengths, in particular of the a^, which is known with less precision. 
These values are, however, subject to further improvement and should be considered 
preliminary. It should also be noted that they are in perfect agreement with the theoretical 
prediction by Caprini et al. of ^JJ^ = 44llg^ — ?'272^jq g. 




The three authors of this work together with F. J. Yndurain (see J. R. Pelaez talk in 
this conference [[8|]) have derived a modified set of Roy-like equations which are based 
on once- subtracted dispersion relations - Roy's equations are twice subtracted. These 
new equations (GKPY for brevity) have a very different dependence on the observables. 
In particular, given the same input, the uncertainty dependence on scattering lengths is 
much weaker. Indeed, at low energy (below ~ 350 MeV) we find the uncertainties for 
GKPY equations bigger than the ones given by Roy's equations, but above ~ 400 MeV 
they are already smaller, as they do not increase with energy. This allows us to obtain 
the position of the sigma pole from Constrained Fits to Data with higher accuracy than 
using standard Roy's equations alone. 

Moreover, we have already performed a preliminary Constrained Fit to Data (CFD-II) 
in which these new equations are also imposed as new constraints within errors (see [8*]). 
This gives rise to a new set of parameters which better encode the experimental infor- 
mation together with unitarity, analyticity and crossing symmetry, therefore allowing us 
to obtain a more precise and reliable determination of the sigma pole. The result for the 
sigma pole position using the very preliminary CFD-II set is: 




(7) 




WORK IN PROGRESS 




(462 ±51)- / (264 ± 20) MeV, (preliminary from Roy Eqs.) (8) 
(458 ± 15) - /(262 ± 15) MeV, (preliminary from GKPY) (9) 



although, as it can be seen in the Figure more work is needed on the /o(980) 
region, as GKPY equations can discern among solutions that were equivalent for Roy's 
equations. The analysis should be complete within the next few months. 




Constrained Fits to Data (FDR+SR+Roy+GKPY) 
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FIGURE 1. a) Difference between input and output of GKPY equations for Re , with the correspond- 
ing error band. The smaller uncertainties carried by these new equations will allow us to improve the fit 
in the /o(980) region, b) Relative difference between the conformal expansion and Roy's equations. The 
contours show the regions inside of which the difference is less than 5 and 10 percent. The sigma pole 
positions with their errors are plotted as ellipses. 



THE CONFORMAL EXPANSION 

The most rigorous way to extrapolate to the complex plane is by using Roy's or GKPY 
equations. However, dealing with the whole set of equations is complicated and compu- 
tationally tedious. For simple applications, there exists a simple approximate solution, 
which is very easy to handle: the conformal expansion. This is a model independent 
parametrization of the experimental data at low energies, based on unitarity and elas- 
ticity, and can describe experimental data with few parameters. The explicit details are 
explained in full length in Ref. [9], here it suffices to remember that, for elastic scat- 
tering, a given partial amplitude of definite isospin / and angular momentum I can be 
written as t'^{s) = with 1/^(5) being the effective range function, which can be 

series expanded in the conformal variable (o(s) = ^7^/^ as follows: 

V^W = ( T^ + Bo + B,CD{^) +B2C0{sf + ..] (10) 

S-Zq/2 [Mny/S J 

Three parameters are enough to describe the experimental data below the nn KK 
inelastic threshold in all studied KK partial waves (5, P, D, F and G). 

We show in Figure [ih the complex plane for ^/s, plotting the contours inside of 
which the relative difference between the amplitude calculated with GKPY equations 

and calculated with the conformal expansion, i.e., Annf*) = ti — 7 t4s less than 

5 or 10 percent. Note that we use our preliminary parameters obtained by constraining 
the data fit to satisfy FDR, Roy's equations and GKPY equations. We can see that for 
the sigma region, the difference between both calculations is less than 5 percent, and 
indeed the pole locations overlap within their errors: using the same input as for Eq. ^ 



we find: 

\J = {A1%±\1)- i {262 ±1) MeV. (11) 

We see that the systematic uncertainties associated with the conformal expansion, which 
correspond to neglecting crossing symmetry and the inelastic cut — /o(980) region — 
are bigger than our previous estimation in Ref. [£], but still within the 5 percent band. 
The crudeness of our estimation of systematic uncertainties has been pointed out by 
I. Caprini in this conference [|llil - Certainly in ref. |@] we only estimated crudely the 
systematic error from different parametrizations and truncations of our conformal ex- 
pansion, and we came up with Ayy^Mfj = ±11 MeV and A.sys^a/'^ = ±2 MeV. I. Caprini 



[|11|1 . using an arbitrary sampling of conformal parametrizations and conformal vari- 
ables - which are different from ours - has provided a new estimate /S.sysM(, ~ 40 MeV 
and A„M,,rc7/2 = ±40 MeV, namely a factor of 3 or 4 larger than ours. Nevertheless, 
the difference between central values for the sigma pole position obtained from GKPY 
equations and our conformal expansion provides the systematic theoretical uncertainties 
in our conformal expansion calculation, without depending on a sampling of conformal 
variables and parametrizations. In view of the differences between Eqs.® and ([TT]) . we 
certainly agree that the systematic uncertainty of the pole obtained from our conformal 
parametrization was underestimated, although not as much as it is suggested in j TH ] but 
just by a factor of two l^sys^a — 20 MeV. 

To summarize this section, up to the sigma region of the complex plane, our conformal 
expansion provides a very simple and reliable approximation, accurate to 5%, to the tzk 
scattering amplitude as obtained from data fits constrained to satisfy FDR, Roy and 
GKPY equations. 
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